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Abstract 

We calculate the electromagnetic form factor of the pion in lattice gauge the- 
ory. The non-perturbatively improved Sheikoleslami-Wohlert lattice action is used 
together with the 0(a) improved current. The form factor is compared to results for 
other choices for the current and features of the structure of the pion deduced from 
the 'Bethe-Salpeter wave function' are discussed. 



PACS Indices: 11.15.Ha, 12.38.Gc 



1 Introduction 



The pion as the simplest particle with only two valence quarks has been the subject of 
many studies. Global features of the pions - their charge and spin - are easily incorporated 
in model calculations. The form factor, which directly reflects the internal structure of this 
elementary particle, is clearly an important challenge. Many earlier calculations are based 
on ad hoc models that model QCD or sum over selected subsets of Feynman diagrams. 
However, the most reliable approach, in particular when addressing non-perturbative fea- 
tures as the electromagnetic form factor at intermediate momentum transfers, is the use 
of lattice QCD. The first lattice results were obtained by Martinelli and Sachrajda pQ, 
which was followed by a more detailed study by Draper et al. 0, who showed that the 
form factor obtained through lattice QCD with the Wilson action could be described by 
a simple monopole form as suggested by vector meson dominance |3]. Below, we extend 
these early studies in two ways. We use an improved lattice action and an 0(a) improved 
electromagnetic current operator. Furthermore, we also extend the calculations to lower 
pion masses than achieved before. Several features of the internal structure of the pion 
have been obtained previously jH 13 IE] by calculating the 'Bethe-Salpeter wave func- 
tion', which can be used to estimate the relative separation of the quark- ant iquark pair 
in the pion. We also use this approach and compare its predictions to the results of our 
direct calculation of the pion form factor. 

2 The method 

In comparison to earlier lattice calculations of pion properties, the major difference of our 
approach is the systematic reduction of the discretisation error in the calculation of the 
matrix elements. We use the non-perturbatively 0(a) improved [Hj clover action JU] and 
the corresponding 0(a) improved current [TTJ [T2J EH] . 

Using this action, we proceed analogous to [2] and calculate the two- and three-point 
correlation functions for the pion. Projecting onto definite pion three-momentum, the 
two-point function is 

G 2 (t,p)=E(^x)0 t (O,O)) e^' x , (1) 

X 

where is the operator projecting on a state with the pion quantum numbers. Below, we 
will consider a 7r + meson, consisting of a u and d quark. Neglecting all spin, colour and 
flavour indices, this operator is given by 

(j)(x) = %p{x) 75 ip(x) . (2) 

In the three point function, which yields the desired form factor, we project onto specific 
initial and final pion three momenta, pj and pj 

G 3 (t f ,t;p f , Pi ) = Yl (0(x / )j4(x)0 t (O))e-P/-( x /- x )-P- x . (3) 

X/ ,X 

This function involves the electromagnetic current operator j^; since here we use only the 
component /x = 4, we do not include a /i-index in the definition of the three point function. 
It is well known that the local current, 

jfa) = $(x) ltl iP(x), (4) 
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is not conserved on the lattice. The conserved Noether current that belongs to our action, 

3° = k (ip(x)(l - ltx )U^x)i){x + fi)- $(x + A)(l + ^)Ul(x)iP{x)) , (5) 

is identical to the conserved current for the Wilson action and still contains corrections of 
O(a) at Q 2 7^ 0; we use Q 2 = —q 2 , where q is the four momentum transfer to the pion. 
The conserved and improved vector current j 1 is of the form |TT|. IT2| IT3*] 

jj, = z v {tf 0) + ac v d l/ T flu } , (6) 

with 

T/jh, = $(x)i (Tfj, v ip(x) , (7) 
Zy = Z v (l + ab v m q ). 

The bare-quark mass is obtained from: 

1/1 1 \ 

am « = 2{-.-^) (8) 

where k c is the kappa value in the chiral limit and a the lattice spacing. For our simulation 
we use k c = 0.13525 [14j . The constants in j 1 are determined such that the matrix element 
of the current operator receives no correction to 0(a). 



3 Details of the calculation 

Our calculations were carried out in the quenched approximation on an 
24 3 x 32 lattice and based on a set of 100 configurations for the link variables at /3 — 6 and 
csw — 1-769 [9j. After an initial thermalisation of 2500 sweeps, we obtained configurations 
at intervals of 500 sweeps. Each sweep consists of a pseudo-heatbath step with FHKP 
updating in the SU(2) subgroups, followed by four over-relaxation steps. In contrast to 
the Dirichlet conditions in [2] , we impose anti-periodic boundary conditions on the quarks 
and periodic boundary conditions on the gluons. Three values of the hopping parameter 
k were used 

«i = 0.1323 , k 2 = 0.1338 , and k 3 = 0.1343. (9) 

For the improved current, we use the parameters Zy, by and cy as determined by Bhat- 
tacharya et al. [T5] . 

Conservation of the total charge generated at the source at t — provides a test |2 3 for 
our calculation, relating the fi = 4 component of the three-point function for q = to the 
two-point function. For our periodic boundary conditions, it reads 

G 3 (t f , t; p, p) - G 3 (t f , t'; p, p) = G 2 {t f , p) , (10) 

where tf < t' < N T . We find that all configurations we use each satisfy this condition to 
at least 1 ppm. 

For the results discussed below, we chose the pion three-momenta such that |p,| 2 = 
\pf\ 2 = 2 in units of the minimal momentum for our lattice. This guarantees for the 
elastic pion form factor that Ef — Ei = q = 0; it greatly simplifies the kinematic factors 
appearing in the three-point function. Different values for the three momentum transfer q 
were obtained by varying the relative orientation of the initial and final pion momenta. 
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In order to improve the projection onto the ground state, we smeared the pion operator 
at the sink in G2 and G3 by the method proposed in [7j. We found that a quark-antiquark 
distance R = 3 works best. The quark-antiquark pair was connected by APE smeared 
gluon links at smearing level 4 and relative weight 2 between straight links and staples. 

To extract the desired information from our numerical results, we assume the two-point 
function of the pion to have the form 

G 2 (t, P) = E y/Z&p) Z "(p) e" * cosh{^ - *)} , (11) 

n=0 Z 

including the contribution of the ground state (n=0) and a first excited one (n=l). The 
Z\ denote the matrix elements, 

Z™(p) = |(fi|0 R |n,p)| 2 , (12) 

and E®, are the energies of ground and excited state, respectively; the subscript R 
indicates the operator smearing. 

The three-point function is parametrised as 

G 3 (t f ,t;p f , Pi ) = F(Q 2 ) JZ%(p f ) Z«{ Vi )e- E *f <«/"*> "'V 
+ {^(P/)^o (pO(l,P/bM(0)|0,P/>e-^/ (t '- t) -^ t +(1-0)} . (13) 

Effects involving, for example, the production of pion pairs, as well as 'wrap around ef- 
fects' due to the propagation of states beyond tf are exponentially suppressed (< C(e~ 5 )); 
similarly, an elastic contribution from the excited state was estimated to be of the order 
of 1% or less. All these effects are not reflected in our chosen parametrisation. 

All parameters in the 2- and 3-point functions - energies E, Z-factors and the form 
factor F(Q 2 ) - were fit simultaneously to the data from all configurations. For the three- 
point function, we chose tf = 11 and let the current insertion time t vary from to 10. For 
maximum spatial symmetry, all values corresponding to the same value |p| in the two-point 
function and all Pij yielding the same q in the three-point function were combined for the 
fit. The value for the parameters and their error in these simultaneous fits was obtained 
through a single elimination jackknife procedure. Since we satisfy eq. (|1U|) to high accuracy, 
we show F(0)=1 in the results below instead of using the result from a fit at Q 2 = 0, which 
would be less accurate in this case. 



4 Results 

Our method to extract the pion form factor is non-perturbatively improved in two respects: 
we use an improved action and an improved current operator. We can get an impression 
of the importance of the latter effect by comparing the conserved Noether current corre- 
sponding to the improved action with the improved current (which is also conserved). The 
results 1 are shown in FigQ for the lightest of our three quark masses. The form factor 
from the improved current is systematically lower than the one from the conserved current. 
The difference grows with Q 2 and reaches about 25% at the largest momentum transfer 
considered here. 

1 In all our results we only show the statistical errors. 
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Figure 1: Form factors for the three currents with lowest quark mass. Solid line: monopole 
form with m p taken from literature (see text). Data for the local current shifted horizontally 
for clarity. 



The structure of the improved current can be further understood by comparing the 
improved current to the renormalised local current, 

ft* = Z V ^ 1 ^, (14) 

which is not conserved. For our kinematics with q = 0, we can also extract a form factor 
from j\. It is also shown in FigQand can be seen to lie very close to the improved current. 
This means that the contribution of the term proportional to cy in j{ is very small. Closer 
inspection shows that while the matrix element of the tensor term can become almost 
comparable to that of the 74 term, the overall tensor contribution is reduced by the small 
coefficient cy. Similar statements also hold for our two additional k- values. 

It is worth mentioning that with the Zy, cy and by values taken from and per- 
forming a fit at Q 2 = we obtained F 1 / F c = 1 to better than 1% with a statistical error 
of about 5%. 

In the previous study [2] of the pion form factor, where the Wilson action was used, 
the results were compared to a monopole form factor 




a form suggested by vector meson dominance. We also show in Fig^ a monopole form 
factor using the value for the p-mass obtained by interpolating the results from ^3] which 
uses the same action as we do. This monopole form factor describes our results for the 
improved current at all but the highest Q 2 very well. As in j2], we observe that the 
conserved current lies consistently above the monopole form factor. A similar behaviour 
was found also for our other two K-values. 

In FigEl we show our results for improved form factors for all three values for k. The 
corresponding quark- and pion masses are given in Tabled The form factors systematically 
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Figure 2: Form factors as a function of Q 2 for the three pion masses. Curves: monopole 
fits to the data. 

decrease with decreasing pion mass. The form factors for the two lightest pions, with 2 
777,^=541 and 671 MeV, come very close together. As can be seen, the statistical error of 
the extracted form factors grows as the quark mass decreases. Nevertheless, we are still 
able to obtain conclusive results for the highest k. The corresponding pion mass of 541 
MeV is substantially lower than in the previous work, where ~ 1 GeV. Given that we 
have only three data sets, we have not attempted to extrapolate our improved form factor 
to the physical pion mass. 

We also fitted our results for the improved form factors to a monopole form factor. 
In doing so, we omitted the highest momentum data point and extracted in each case a 
vector meson mass, my, shown in Tabled They are close to the values for m p taken from 
interpolations to literature data [Tlj . 

In examining the two-point Green function for various quark-antiquark distances, we 
also obtain the 'Bethe-Salpeter wavef unction', 



RK ' (16) 



Following the procedure in 0IB], we obtain (r 2 )g S , shown in Tabled These RMS-radii 
are compared to the values extracted from the low-Q 2 behaviour of the form factor, 



alQ 2 



l -(r 2 ) =-\ (17) 
6 \ iff ml y ' 



where in the last step we have assumed a monopole form and use the fitted parameter 
my. In agreement with the findings of [UIHIIE], we see that the Bethe-Salpeter predictions 

2 For definiteness, we have taken the lattice spacing a — 0.105 fm from 
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0.13230 
0.13380 
0.13430 


0.082 
0.040 
0.026 


0.515(2) 
0.357(2) 
0.288(2) 


0.625(5) 
0.513(5) 
0.476(7) 


0.597(14) 
0.496(15) 
0.470(19) 


2.530(2) 
2.615(2) 
2.629(2) 


4.23(10) 
4.94(15) 
5.21(21) 



Table 1: Masses and RMS-values for the different kappa values, in lattice units. The 
p-mass has been taken from \14j - 

are very insensitive to the value of the quark or pion mass. However, it is well known |Hj 
that the information that can be obtained from the Bethe-Salpeter approach as described 
above is only an approximation. It assumes, in the extraction of (r 2 ), that the center of 
mass of the pion is always halfway between the valence quark and antiquark, not allowing 
for the motion of the gluons. The extraction of (r 2 ) from the calculated pion form factor 
does not involve this restriction for the valence (anti-) quark motion. As can be seen, the 
more reliable determination from F(Q 2 ) leads, as expected, to a larger radius. Moreover, 
this radius shows a substantial dependence on the mass. 

We have presented here the first calculation of the electromagnetic form factor of the 
pion based on an 0(a) improved action and the concomitant improved vector current. This 
is seen to lead to significant changes in the prediction for the internal structure of the pion. 
We observe a decrease of the form factor for decreasing pion mass, which in turn leads 
to an increase of the RMS-radius. This mass- dependence of the radius is not seen in the 
Bethe-Salpeter approach. Furthermore, the mass of the pion we reach in our calculations 
is considerably closer to the physical value than in previous work. 

The computational effort involved in taking the improvement into account is small. 
Since it guarantees elimination of 0(a) discretisation errors to all orders in the coupling 
constant, use of this method in future work seems logical in pushing the calculations further 
towards the physical limit. 
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